Quantum algorithm for finding periodicities in the spectrum 
of a black-box Hamiltonian or unitary transformation 
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Estimating the eigenvalues of a unitary transformation U by standard phase estimation requires 
the implementation of controlled-C/-gates which are not available if U is only given as a black box. 

We show that a simple trick allows to measure eigenvalues of U ® If' 1 even in this case. Running 
the algorithm several times allows therefore to estimate the autocorrelation function of the density 
of eigenstates of U. This can be applied to find periodicities in the energy spectrum of a quantum 
system with unknown Hamiltonian if it can be coupled to a quantum computer. 



I. STANDARD PHASE ESTIMATION AND ITS 
WEAKNESS 

Finding the eigenvalues of unitary transformations or 
self-adjoint operators is a central task in quantum me- 
chanics. The thermodynamic and dynamical properties 
of a quantum system are determined by the spectrum of 
a Hamiltonian and the corresponding unitary transfor- 
mations. Furthermore the estimation of eigenvalues is 
an important tool in quantum computation (see |IJ). For 
that reason the algorithm for phase estimation has been 
developed (see and references therein). We rephrase 
it as follows. We have a Hilbert space lZ a <g> Ti where Ti 
is the target register where the considered unitary U acts 
on and an ancilla register lZ a consisting of k qubits if an 
accuracy of the eigenvalues of U of the order 2~ k is de- 
sired. Assume the target register to be in an eigenstate 
\^pu) °f U with eigenvalue exjp(icf>). Initialize the ancilla 
register in an equal superposition of all its logical states, 
i.e., 
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where I is the binary number corresponding to the k an- 
cilla qubits j = 0, . . . , k — 1. On the joint Hilbert space 
Ha <£> Ti apply for all j = 0, . . . , k — 1 the transformations 



lj)(lj\®U 23 + \0j)(0j\®l, 



(1) 



where and |0j)(0j| are the projectors onto the 

|0) and |1) state of the ancilla qubit j, respectively. The 
operation U 2 ' is the 2 J -fold iteration of U. 

If TL is in the state \tpn) then \ijjn) is converted into 
the state 

by the 'kick-back-effect' Q. After Fourier transformation 
on lZ a we obtain 
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i.e., the probability distribution is peaked around m = 
4>2 k /ir. If \ipn) is not an eigenstate of U, than the algo- 
rithm will project approximatively (in the limit of large 
k) onto any of the eigenstates ||. If the initial state on 
TL is a density matrix which is diagonal in the basis of 
U, one will obtain any of the eigenvalues of U with the 
corresponding probability. 

At first sight, quantum phase estimation seems to be 
applicable for finding energy values and eigenstates of an 
unknown Hamiltonian of an arbitrary quantum system 
simply by setting U := exp(—iHt). This would be inter- 
esting for the investigation of complex physical systems. 
To measure eigenvalues of interaction Hamiltonians in 
many-spin systems, as molecules or solid states, for in- 
stance, would be rather useful. 

But there is a severe problem which is essentially that 
quantum phase estimation does not use the implementa- 
tion of U as a black box subroutine. It uses the conditional 
transformations Vj (see eq. (Q)) as black boxes and one 
should emphasize that no canonical conversion procedure 
building Vj from U is known if U is a black box. Some- 
times this fact is hidden by using a language too classical 
if one explains the action of Vj by claiming that it im- 
plements U (or its iterations) if the corresponding ancilla 
qubit is in the logical state |1). This hides the fact that a 
superposition state of the ancilla has to lead to a super- 
position of the two actions 'implementation of IP and 'no 
implementation'. Imagine that the black box implement- 
ing U contains a memory storing the information whether 
U is implemented or not. Then such a superposition of 
implementation is destroyed by the memory. However, 
in Q| we have shown that the unitary evolution U can 
in principle be conjugated by other unitary transforma- 
tions in such a way that the net effect is a ' controlled- IP. 
But the class of unitary transformations U considered 
there is only the evolutions according to rt-qubit pair- 
interaction Hamiltonians. Here we address the question 
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how to use quantum phase estimation for obtaining in- 
formation about the spectrum of U if we have no prior 
information about U at all. We show that it is at least 
possible to get the autocorrelation function of the spec- 
trum of U, or, speaking more explicitly, the spectrum of 
U<8)U\ provided that the following assumptions are true: 

1. The operation U is implementable on a system 
Ti. which can be brought into interaction with an- 
other register TZi of equal Hilbert space dimension 
in such a way that complete exchange of quan- 
tum information between TL and TZi is possible. 
If U = exp(— iHt) for an appropriate t > and 
H is the real Hamiltonian of the system TL we as- 
sume that this exchange of information can either 
be done on a small time scale compared to the evo- 
lution according to H or the natural evolution can 
be switched off during the implementation of this 
information exchange. 

2. There is another quantum register TZ2 with the 
same dimension as TL and TZi and an ancilla regis- 
ter TZ a consisting of k bits if the desired accuracy 
for the eigenvalues is 2 k . 

3. On the system TZi ®7Z2®lZ a we have a set of quan- 
tum transformations available which is universal for 
quantum computation. 

Of course the assumption that the unknown Hamil- 
tonian H can be switched off is problematic, but if ad- 
ditional prior information about the structure of H is 
available, standard decoupling techniques || can be used. 
Note that assumption 1 is considerably weaker than the 
assumption that H can be switched on and off by the 
quantum state of an ancilla qubit]^. 

II. IMPLEMENTING A CONDITIONAL 
TRANSFORMATION 

The essential part of our algorithm is rather simple in 
contrast to Q for the cost that we obtain eigenvectors 
and eigenvalues of U (g> U~ x instead of those of U. It con- 
sists of a conjugation of U by known unitary transforma- 
tions in such a way that the net effect is the conditional 
transformation 

V! := (8 U v <g> 1 + IOjXOjI <8» 1 O U v , 

where |0j) and are states of the ancilla qubit j. One 
can see easily that the effect on the ancilla's states is the 
same if the unconditional unitary U~ 21 is implemented 
on TZ2 after each implementation of Vj. This implements 
the conditional transformation 



J This is discussed in more detail in 



V? := ® U 2J ® U- 2 ' + \0 3 )(0 3 \ ® 1 ® 1. 

Using standard phase estimation we can use V" for 
obtaining eigenvalues of U ®U*. 

The procedure for implementing V- consists of the fol- 
lowing steps for j = 0, . . . , k — 1. 

1. Implement state exchange of the registers TL and 
Til, i.e., the unitary W with W(\a) ® \(3)) := 

®|a). 

2. Implement U 2 ' on TL. If U = exp(-iHt), i.e., if 
u is the unitary evolution according to the system 
Hamiltonian H, then one has to wait the time 2H. 

3. Implement state exchange of TL and TZ\ again. 
Steps 1-3 implement the transformation exp(-iHt) 
on the register TZi. 

4. Implement a conditional exchange of the states of 
TZi and TZ2 depending on the state of qubit j in the 
ancilla state, i.e. we implement 

\lj)(lj\®W + 10^18 1 

on TZ a ®TZ\® TZi- 

The conditional exchange can easily be implemented, 
if the registers TZ\ and TZi consist of qubits. In this case 
the conditioned permutation of two corresponding qubits 
is a usual Fredkin-gate Q]. 

Our algorithm might have applications for investigat- 
ing the spectrum of a many-particle Hamiltonian in solid- 
state physics, since the spectrum and its gaps deter- 
mines dynamical and thermo-dynamical behavior of the 
system ||. Of course for many-particle systems it is 
not possible to find the complete set of eigenvalues of 
cxp(-iHt) (g) exp(iHt) since the dimension of Ti grows 
exponentially. But periodicities in the spectrum of H 
could be found. We sketch this idea. First we assume 
that t is chosen in such a way that tA < n, where A 
is an upper bound on the difference between the largest 
and smallest eigenvalue of H given by prior knowledge. 
Assume that TZi and 7Z2 are prepared in the same initial 
density matrix p. The state p defines a probability mea- 
sure on R by setting p(A) = if A is no eigenvalue of u and 
p(X) := (X\p\X) else. If p is the maximally mixed state 
and the dimension of TL is large such that the probability 
measure can approximatively described by a probability 
density, this measure is known as the density of states |Q . 
Running the algorithm several times allows to estimate 
the density of eigenstates of H® 1 — l<g>H which is the au- 
tocorrelation function of the density of eigenstates of H . 
If this density contains periodicities, for instance when 
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there are spectral gaps with equal distances, they can be pend on n. Hence the size of the required ancilla register 
detected by using our algorithm. Note that in solid-state does not necessarily grow with n for detecting interesting 
physics, for instance, energy gaps occur which do not de- gaps. 
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